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We have studied the nonhnear dynamics of whistler waves in magnetized plasmas. Since plas- 
mas and beam-plasma systems considered here are assumed to be weakly collisional, the point of 
reference for the analysis performed in the present paper is the system of hydrodynamic and field 
equations. We have applied the renormalization group method to obtain dynamical equations for 
the slowly varying amplitudes of whistler waves. Further, it has been shown that the amplitudes 
of eigenmodes satisfy an infinite system of coupled nonlinear Schrodinger equations. In this sense, 
the whistler eigenmodes form a sort of a gas of interacting quasiparticles, while the slowly varying 
amplitudes can be considered as dynamical variables heralding the relevant information about the 
system. An important feature of our description is that whistler waves do not perturb the initial 
uniform density of plasma electrons. The plasma response to the induced whistler waves consists in 
velocity redistribution which follows exactly the behaviour of the whistlers. In addition, selection 
rules governing the nonlinear mode coupling have been derived, which represent another interesting 
peculiarity of our description. 
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I. INTRODUCTION 

Over four decades passed since it was first shown that 
plasmas and beam-plasma systems immersed in an ex- 
ternal magnetic field can support travelling electromag- 
netic waves with specific features. These waves propa- 
gate parallel to the applied magnetic field being circu- 
larly polarized in a plane transverse to the direction of 
propagation. It has become conventional in the physics 
of magnetized plasmas to call such structures waves in 
the whistler mode. 

Although the linear stability properties of the electro- 
magnetic waves in the whistler mode are relatively well 
studied P, 0, 0, Hi there is a serious gap in the un- 
derstanding of their nonlinear behaviour. Chen et al. 

have shown that electromagnetic whistler waves can 
be considered as complementary to the nonlinear travel- 
ling electrostatic waves, known as the Bernstein-Greene- 
Kruskal (BGK) modes 0. While the BGK modes are 
longitudinal, the whistler modes are transverse, in other 
words, the components of the electric and magnetic field 
of the whistler wave parallel to the external magnetic field 
are both zero. The study of the nonlinear behaviour of 
whistler waves has been initiated by Taniuti and Washimi 
0, who obtained a nonlinear Schrodinger equation for 
the slowly varying amplitude (see also Reference Q). 
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The present paper is aimed at filling the gap in the un- 
derstanding of the nonlinear evolution of whistler waves. 
The method adopted here is the renormalization group 
(RG) method [3, [l3 ■ The basic feature of this approach 
is that it provides a convenient and straightforward tool 
to obtain an adequate description of the physically es- 
sential properties of self-organization and formation of 
patterns in complex systems. Coherent structures which 
result from the nonlinear interaction between plane waves 
evolve on time and/or spatial scales comparatively large 
compared to those the fast oscillations occur. The RG 
method can be considered as a powerful systematic pro- 
cedure to separate the relatively slow dynamics from the 
fast one, which is of no considerable physical relevance. 
In a context similar to that of the present paper, it has 
been successfully applied by one of the authors [lOlllll l to 
study collective effects in intense charged-particle beams. 

The paper is organized as follows. In the next section, 
we state the basic equations which will be the subject of 
the renormalization group reduction in section III. Start- 
ing from a single equation [see equation ljIII.6|) ] for the 
electromagnetic vector potential, we obtain a formal per- 
turbation expansion of its solution to second order. As 
expected, it contains secular terms proportional to pow- 
ers of the time variable which is the only renormaliza- 
tion parameter adopted in our approach. In section IV, 
the arbitrary constant amplitudes of the perturbation ex- 
pansion are renormalized such as to eliminate the secular 
terms. As a result, a set of equations for the renormal- 
ized slowly varying amplitudes is obtained, known as the 
renormalization group equations (RGEs). These equa- 
tions comprise an infinite system of coupled nonlinear 
Schrodinger equations. In section V, the latter are ana- 
lyzed in the simplest case. Finally, section VI is dedicated 
to discussion and conclusions. 
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II. FORMULATION OF THE PROBLEM AND 
BASIC EQUATIONS 

Plasmas and beam-plasma systems considered in the 
present paper are assumed to be weakly coUisional. 
Therefore, the dynamics of plasma species is well de- 
scribed by the hydrodynamic equations coupled with the 
equations for the electromagnetic self- fields. We start 
with the equations for plasma in an external constant 
magnetic field Bq, which can be written as follows 



^ + V • {naVa) = 0, 



(III) 



where e is a formal small parameter introduced for con- 
venience, which will be set equal to one at the end of the 
calculations. Thus, the basic equations to be used for the 
subsequent analysis can be written in the form 



^ + noV • V + eV • (iVV) = 0, 



^ + .V.VV^-^M_viV 

dt m{7io + eN) 



(11.10) 



DgVa 



maUa 



V, X (Bo + B)], 



(II.2) 

where and are the density and the current veloc- 
ity of the species a. Furthermore, rua, qa and Ta are the 
mass, the relative charge and the temperature, respec- 
tively, while fcs is the Boltzmann constant. The substan- 
tional derivative on the left-hand-side of equation l|II.2p 
is defined as 



(II.3) 



The electromagnetic self-fields E and B can be obtained 
in terms of the electromagnetic vector A and scalar if 
potentials according to the well-known relations 



dA 



E=-V(/?-^, B = V X A. (II.4) 
The latter satisfy the wave equations 



□ A = -Moe ^ Wa^a Va , 0(p 



e \ ^ 



£0 



naQa, 



in the Lorcntz gauge 



(II.5) 



(II.6) 



Here □ denotes the well-known d'Alembert operator. 
In what follows, we consider the case of a quasineutral 
plasma 



naQa = 0, 



(II.7) 



in a constant external magnetic field along the cc-axis 
Bq = (i?o,0, 0). Then, equations HII.l|) - (|II.6p possess a 
stationary solution 

n-a = fiaQ = const, Va = 0, A = 0, (p = 0. (II. 8) 

The frequency of the wave will be taken as much higher 
than the ion-cyclotron frequency. Therefore, we can fur- 
ther neglect the ion motion and scale the hydrodynamic 
and field variables as 



-[E + Vx (Bo-^eB)], 



(11.11) 



□A = /ioe(no + eiV)V, ^^ + V-A = 0. (11.12) 



Before we continue with the renormalization group re- 
duction of the system of equations l|II.10(l - (|II.12(l in the 
next section, let us assume that the actual dependence 
of the quantities N, V, A and ip on the spatial variables 
is represented by the expression 



*-*(x,X;i), * = (7V,V,A,^), 



(11.13) 



where X = ex is a slow spatial variable. Thus, the only 
renormalization parameter left at our disposal is the time 
t which will prove extremely convenient and simplify te- 
dious algebra in the sequel. 



III. RENORMALIZATION GROUP 
REDUCTION OF THE 
MAGNETOHYDRODYNAMIC EQUATIONS 

Following the standard procedure of the renormaliza- 
tion group method, we represent ^' as a perturbation 
expansion 



n=0 



(IILl) 



in the formal small parameter e. The next step con- 
sists in expanding the system of hydrodynamic and field 
equations (|II.10|l - ljII.12|l in the small parameter e, and 
obtaining their naive perturbation solution order by or- 
der. Note that in all orders the perturbation equations 
acquire the general form 



dt 



(IIL2) 



no + eiV, Ve = eV, 



eA, 



(11.9) 



dt 



no 



-E„ 



UcV„ X e, + W„, (IIL3) 
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□A„ = /ioenoV„ + U„, 



1 dipn 



dt 



(III.4) 

where a„, /3„, U„ and W„ are quantities, that have been 
already determined from previous orders. Here 



kaT eBn 

uJc=— III.5 
m m 



are the thermal velocity of electrons and the electron- 
cyclotron frequency, respectively and Bx = (1,0,0) is the 
unit vector in the z-direction. Manipulating in an ob- 
it is possible to 



vious manner equations ljIII.2|) - (|III.4l 
obtain a single equation for A„. The latter reads as 



□ 



■ 4nv(v • A„) 



n^A,, 
at 



- WpV(V • A„) = ^oeno 



9W„ 

dt 



-^oev^Van — V(V ■ U„) 



dVn 



dt 



9*2 



where 



(III.6) 



(III.7) 



is the electron plasma frequency. Note that the thermal 
velocity vj- as defined by equation H111.5() can be alterna- 
tively expressed according to the expression 



(III., 



where rjj is the electron Debyc radius. Equation lllll.6|) 
represents the starting point for the renormalization 
group reduction, the final goal of which is to obtain a 
description of the relatively slow dynamics leading to for- 
mation of patterns and coherent structures. 

Let us proceed order by order. We assume that the 
dependence on the fast spatial variables x = (x, y, z) is 
through the longitudinal (parallel to the external mag- 
netic field Bp) x-coordinate only. The solution to the 
zero-order perturbation equations ljIII.6|l can be written 
as 



where 



Ao=^A(°)Ae' 



il^k{x] t) ^ kx - ujkt, 



(III.9) 



(III.IO) 



and Ak is an infinite set of constant complex amplitudes, 
which will be the subject of the renormalization proce- 
dure in the sequel. Here "constant" means that the am- 
plitudes Ak do not depend on the fast spatial variable x 



and on the time t, however, it can depend on the slow spa- 
tial variables X. The summation sign in equation (|111.9() 
and throughout the paper implies summation over the 
wave number k in the case where it takes discrete values, 
or integration in the continuous case. From the disper- 
sion equation 



V{k\LOk) = Lol 



0, (HI.ll) 



it follows that the wave frequency uJk can be expressed 
in terms of the wave number fc, where the Fourier-image 
□fc of the d'Alembert operator can be written according 
to 



(HI.12) 



Moreover, it can be verified in a straightforward manner 
that the constant vector A^"-* can be expressed as 



A^°^ = (0,l,-*sgn(fc)), 



(IH.13) 



where sgn(fc) is the well-known sign-function. Details 
concerning the derivation of the dispersion law Ijlll.lip 
and equation (jllLlSp can be found in the Appendix. Note 
that equation ljIII.13|) is an alternative representation of 
the solvability condition (|A.3|I . It is important to em- 
phasize that 



tO-k 



A. 



A 



(HI. 14) 



where the asterisk denotes complex conjugation. The lat- 
ter assures that the vector potential as defined by equa- 
tion ljIII.9p is a real quantity. The zero-order current ve- 
locity Vq obtained directly from the first equation (|III.4I) 
can be written as 



Vo 



V 



(0) 



(IH.15) 



In addition, the zero-order density, scalar potential and 
magnetic field are represented by the expressions 



No = 0, 



^0 = 0, Bo = ^Bi°Ufce''^^ (IH.16) 



where 

B^"'> = -fcA^°^sgn(fc) = (0, -ksgnik),ik). (IH.17) 

It has been mentioned that the first-order "source 
terms" on the right-hand-side of equation (|III.6|I can be 
expressed via quantities already known from zero order. 
Thus, we have 



Ui = -2V • VAo, Wi 



-V-An 



(IH.18) 



m 



-Vo X Bo, (IH.19) 
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where the shorthand notation 



_d_ 



(III.20) 



has been introduced. Note that the vector Wi repre- 
senting the zero-order Lorentz force has the only nonzero 
component along the external magnetic field, that is 



fc.i 



(III.21) 



where 



aki 



(kOi + ink)[l - sgn(fc)sgn(/)]. (III.22) 



Equation ljIII.6|l has now two types of solutions. The 
first is a secular solution linearly dependent on the time 
variable in the first-order approximation. As a rule, the 
highest power in the renormalization parameter of the 
secular terms contained in the standard perturbation ex- 
pansion is equal to the corresponding order in the small 
perturbation parameter. The second solution of equa- 
tion (|III.6|I arising from the nonlinear interaction between 
waves in the first order, is regular. Omitting tedious but 
standard algebra, we present here only the result 



k,l 



where 



Ai^^ = (A«,^^L^-*i^l>gn(fc)j, 



r(i) 



r(i)c 



(III.23) 
(III.24) 



Some of the details of the calculations are presented in 
the Appendix. In explicit form, the components of the 
vector operator A^^"* and those of the infinite matrix A^.;"* 
are given by the expressions 

ik(3k 



,(0) 



(III.25) 



In addition, the constants ak, (3k, Ik and Xk entering the 
expressions above are given by 



o-k = Dfe 



? 9 



Xk = ^k + 



' D 



(III.31) 



(III.32) 



Furthermore, the first-order current velocity can be ex- 
pressed as 

= 51 vi'Ufee*'^^ + 51 AkAi&'^^-^'^'\ 

k k,l 

(III.33) 

where 

Vi'^ = {V^'J , V^'J , ~^V^'Jsgnik)) . (III.34) 

The corresponding operators and matrix coefficients can 
be written explicitly according to the expressions 



(1) 



(1) 



Moeno 



1 



^kx ' ^kl 



(1) 



□ 



kl ,(1) 



A):>, (III.35) 



ky 



X 



(III.36) 

Calculating the first-order density iVi from equation 



(IIII.2II. we obtain 



iii'k+'tf'l) 



N, 



(1) 



□ 



Hoeojk 



ujk V '"^ 



iVk 



(III.37) 



(III.38) 



^"^^ - -2u;kakSgn{k)+^cXk' <^ = z^o^^P.. + " ^g"(')^Sn(0]- 

(III.39) 

ji) e ujk+i^i,,„ ,„ Analogously, for the first-order scalar potential cpi, we 

= 2^d;:;v^^^°' + - find 

(HI 27) 

where ' =Y.^i''> Ake^"^" +Y,vi';> AkAie'^''''+'^'\ (III.40) 



Fk = 2fcwfc[wfcSgn(fc) -I- wjVx, (III.28) 



Uki= ^ ' -{k + lf, (III.29) 



k,l 



(1) - _^7V^^^ - — 



(fc4i'-*V.), (III.41) 



^ 1;^^^+^ _ (fc + ;)2 _ 1 . (ni.30) 2mi;^n,;I?,; 



'Z3 



III.42) 
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Finally, the first-order magnetic field is calculated to be 
Bi =^Bi^)Ae^^^ (III.43) 

k 

where 

B« = (-zsgn(A:)Vfc, -zBf^)sgn(fc)) , (III.44) 



The latter imposes additional restrictions on the depen- 
dence of the wave amplitudes Ak on the slow transverse 
independent variables Y and Z. Straightforward calcu- 
lations yield (see the Appendix for details) 



isgn(fc) 



^(2) ^ 



ky 



(IV.5) 



B 



(1) 

ky 



-Sgn(fc)( 



tkA 



(1) 

ky 



iV 



X 



(III.45) 



A couple of interesting features of the zero and first- 
order perturbation solution are noteworthy to be com- 
mented at this point. First of all, the zero-order density 
Nq vanishes which means that no density waves are in- 
duced by the whistler eigenmodes. The second terms in 
the expressions for the first-order density A^i and current 
velocity Vi [see equations (|III.33ll and HIII.37P ] imply 
contribution from nonlinear interaction between waves 
according to the nonlinear Lorentz force. It will be shown 
in the remainder that these terms give rise to nonlinear 
terms in the renormalization group equation and describe 
solitary wave behaviour of the whistler mode. 



IV. THE RENORMALIZATION GROUP 
EQUATION 

Passing over to the final stage of our renormalization 
group procedure, we note that in second order the quanti- 
ties U2 and W2 entering the right-hand-side of equation 
(IIII.6II can be written as 



-2V • VAi - V^Ao + /ioeTViVo, 



(IV.l) 



W2 = — V(^i-^V7Vi-Vi-VVo-— VixBo, (IV.2) 

m TiQ m 

Since we are interested only in the secular terms in sec- 
ond order, appearing in the expressions for the y and z 
components of the electromagnetic vector potential A2, 
contributions in the source vectors U2 and W2 leading 
to such terms are sufficient for completing the renormal- 
ization group procedure. Thus, we can write 



E 



tA 



(2) 



tJ^'Dr-^ke''" +tT.^ki\Ai\^Ake 



iipk 



(IV.3) 



k,l 



An important remark is in order at this point. From the 
solvability condition (jA.3|l it follows that the complex 
amplitude Ak must satisfy the complex Poisson equation 



(IV.4) 



1 



k^Pk 



^(2) ^ wj.^fcSgn(fc) 
''^ 2wfeafeSgn(fc) 4- ujcXk \ ' Ik^k 



Cky — ~A 



17(1)2 
2^ky . 



where 



,(2) 



OLk 



sgn(fc). 



Gk = Ci^/cSgn(fc)[wfcSgn(fc) -I- Wc]V^ 



VyVfc, 

(IV.6) 
(IV.7) 

(IV.8) 
(IV.9) 



The matrix coefficient Tkiy determining the nonlinear 
contribution represented by the second term in equation 
(|IV.3p reads explicitly as 

1 - sgn(A;)sgn(;) iuJkUi{kUi + lUk)s.gn(k) 



Tk. 



kly 



fionomv^ujiVki 2tjfeafcSgn(fc) + cUcXk 



X [uic{luJk - kuji)sgn{l) + {k + l)iLJkiLJi] 



(IV.IO) 



Following the standard procedure [lyl of the RG method, 
we finally obtain the desired RG equation 



dAk 
dt 



A 



(2) 
ky 



ky 



D):!]Ak 



kly 



Ai 



Ak, (IV. 11) 



where now Ak is the renormalized complex amplitude 
[To|. Thus, the renormalized solution for the electromag- 
netic vector potential acquires the form 



A = 5:A("'Ae 



(IV. 12) 



Analogously, for the electric and magnetic field of the 
whistler wave, one can obtain in a straightforward man- 
ner the following expressions 



B = Y^ B^^^Ake^^^ , E = z ^ cokA^^^Ake^^^ . 

k k 

(IV. 13) 

It is important to mention that the plasma density 
remains unchanged (TV = 0) contrary to the case of elec- 
trostatic waves, where the evolution of the induced elec- 
trostatic waves follows the evolution of the density waves. 
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V. SYSTEM OF COUPLED NONLINEAR 
SCHRODINGER EQUATIONS 

The simplest case of the vahdity of the solvabihty con- 
dition ljIV.4|) consists in the assumption that the slow 
wave amplitudes Ak do not depend on the transverse co- 
ordinates. Setting e = 1 in equation IjlV.lip . we obtain 
the following system of coupled nonlinear Schrodinger 
equations 



isgn(fc) 



dAk 

dt ' 



'WfcSgn(fc) 



dAk 
dx 



A, 



where for simplicity the tilde-sign over the renormalized 
amplitude has been dropped. Moreover, the coefficients 
Vk, Afc and ^ki are given by the expressions 



2A:u;fe[wfcSgn(A:) ■ 



- UJc 



2LJkakSgn{k) + uJcXk ' 



a;fc[cjfcSgn(fc) + ujc] 
2uJkakSgn{k) + uj^Xk 



(V.2) 



4fc^cjfc/3f ^ [cjfcSgn(fc) + UJc 
[2cjfeafeSgn(fc) -I- LUcXkf 



sgn(fc) , (V.3) 



IJ-kl 



1 - sgn(fc)sgn(/) LOkUiikUi + mk) 
fj.onQmv^uJiT>ki 2uJkakSgTi{k) + uj^Xk 



X [ujc{lujk — kuji)s.g\i{l) + (fc + l)ujkijJi\ 



(V.4) 



Interestingly enough, the infinite matrix of coupling co- 
efficients /ife; represents a sort of selection rules. Clearly, 



^J'kk — 0, 



^^k,-k — 0, 



(V.5) 



and 



Hki = 0, for sgn(A:)sgn(0 = 1. (V.6) 



This means that a generic mode with a wave number 
k cannot couple with itself, neither can it couple with 
another mode with a wave number of the same sign. Note 
that this feature is a consequence of the vector character 
of the nonlinear coupling between modes and is due to 
the nonlinear Lorentz force. Therefore, for a given mode 
k the simplest nontrivial reduction of the infinite system 
of coupled nonlinear Schrodinger equations consists of 
minimum two coupled equations. 

Without loss of generality, we can assume in what fol- 
lows that the sign of an arbitrary mode k under consid- 
eration is positive (fc > 0). Suppose that for a particular 
whistler mode with a positive wave number k there exist 



a mode with wave number — / for which the coupling co- 
efficient Hk,-i is maximum. Neglecting all other modes 
but the modes k and —I, we can write 



.dAk , . dAk . d^Ak 



dt 



IVk- 



dx 



dx"^ 



+ liMiV^k, (V.7) 



dt dx dx"^ 



-ti2\Ak\ Ai 



(V.8) 



'Y,^^kl\Al\'Ak, where 
(V.l) 



Ml 



LUkUi{kUi - lUk) 



fj,onomv^u!iVk-i 2ujkak + oJcXk 



X [uJc{kLUl - lujk) + (fc - l)uJk^^i] 



(V.9) 



M2 



iUiUkikUi ~ lUk) 



fionomv^LJkVk-i 2uJiai + uJcXi 



x[LUcikLUi - lujk) + (fc - l)uJku;i] 



(V.IO) 



The system of coupled nonlinear Schrodinger equations 
(|V.7p and IjV.Sp is non integrable in general It rep- 
resents an important starting point for further investiga- 
tions on the nonlinear dynamics and evolution of whistler 
waves in magnetized plasmas. 



VI. DISCUSSION AND CONCLUSIONS 

We studied the nonlinear dynamics of whistler waves 
in magnetized plasmas. Since plasmas and beam-plasma 
systems considered here are assumed to be weakly colli- 
sional, the point of reference for the analysis performed in 
the present paper is the system of hydrodynamic and field 
equations. We apply the renormalization group method 
to obtain dynamical equations for the slowly varying am- 
plitudes of whistler waves. As a result of the investiga- 
tion performed, it has been shown that the amplitudes 
of eigenmodes satisfy an infinite system of coupled non- 
linear Schrodinger equations. In this sense, the whistler 
eigenmodes form a sort of a gas of interacting quasiparti- 
cles, while the slowly varying amplitudes can be consid- 
ered as dynamical variables heralding the relevant infor- 
mation about the system. 

An important feature of our description is that whistler 
waves do not perturb the initial uniform density of 
plasma electrons. The plasma response to the induced 
whistler waves consists in velocity redistribution which 
follows exactly the behaviour of the whistlers. Another 
interesting peculiarity are the selection rules governing 
the nonlinear mode coupling. According to these rules 
modes with the same sign do not couple, which is a direct 
consequence of the vector character of the interaction. 
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Careful inspection shows that the initial source of the 
nonlinear interaction between waves in the whistler mode 
is the zero-order Lorentz force [see equation 
Since the quantity Wi is proportional to A^' 
the above mentioned selection rules follow directly, pro- 
vided the only case in which the cross product does not 
vanish is the case, where modes k and I have different 
sign. 

We believe that the results obtained in the present 
paper might have a wide class of possible applications 
ranging from laboratory experiments to observations of 
a variety of effects relevant to space plasmas. 
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APPENDIX A: DETAILS CONCERNING THE 
DERIVATION OF THE PERTURBATION 
EXPANSION 

Under the assumption that whistler waves propagate 
parallel to the external magnetic field Bq (dependence on 
the longitudinal x-coordinate and the time t only), the 
equation ljIII.61) for the longitudinal component Anx de- 
couple from the equations for the other two components. 
In zero order it has an obvious solution Aqx = 0. The 
perturbation equations for the transverse components of 
the electromagnetic vector potential can be written as 



of equation ljIII.6|) . These yield a secular solution to the 
perturbation equations linear (in first order) with respect 
to the time t. Taking into account the resonant terms 
alone, it can be verified that the operators v4^^^ and 
satisfy the following system of linear equations 



2iu!kakA 



(1) 

ky 



-isgn{k)Fk, 



2iuJkakA^kz 



(A.4) 



(A.5) 



The solvability condition ljA.3|l which also holds in first 
order yields immediately equations (|III.24p and (|III.26p . 
Non resonant terms can be handled in a straightforward 
manner, yielding the second term on the right-hand-side 
of equation ljIII.23|) . 



In second order the right-hand-side of equation (|III.i 
contains resonant terms proportional to e^'^'' , as well as 
terms proportional to te^'^'' . Contributions to the second- 
order solution of the first type can be handled in a way 
similar to that already discussed in first order. To deal 
with the second type of resonant contributions, we write 
the second-order perturbation equations as 



□ 



d^A2y 



dA 



22 



c^; d^A 



2y 



at 



^^sgn(fc)i^fcAe''^^ 
(A.6) 



□ ^r-T^ WcD 



9t2 



^^^^ ^"^'"^'^-tHkAke^K (A.7) 



dt 



-2 



□ 



dt 



to; 92 A 



Oy 



It is straightforward to verify that the solution to equa- 
(^■1) tions l|X6|) and (fO)l is of the form 



□ 



5^ 



OA. 



Oy 



dt 



4 d^Ao, 
'? dt^ 



(A.2) 



To solve the above system of equations, we use the ansatz 
(|III.9|I . Then, the dispersion equation Ijlll.lll) can be 
obtained in a straightforward manner as a condition for 
vanishing of the determinant of the linear system consist- 



ing of the components A^j^^ and A'^^J . Moreover, from the 
dispersion equation it follows that A^^ is proportional to 
that is 



(0) 



A2y,^ = {t^C'ky^. + tLky,,)Ake'^\ (A.8) 
where the coefficients Cky,z and L^y^z can be written as 



Cky — 



2(2wfcQ;fcSgn(fc) + uJcXk) 



■'ky 



. 2/3f)Cfc^sgn(fc) 
'2[jfcafcSgn(fc) -I- LUcXk ' 



(A.9) 



(A.IO) 



A 



(0) 



-*sgn(fc)Ag. 



(A.3) 



Equation (|A.3|I should be regarded as a solvability con- 
dition and must be satisfied order by order. 

The first order is characterized by the presence of reso- 
nant terms (proportional to e'^'^'' ) on the right-hand-side 



Ckz = -iCkySgn{k), Lkz = -iLkySgi\{k). (A. 11) 

Substituting the appropriate form of and collecting 
similar terms proportional to t and , we readily obtain 
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